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We analyze the influence of different bilinear terms giving rise to time reversal invariant topological 
insulating phases in mono and bilayer graphene. We make use of the effective action formalism to 
determine the dependence of the Chern Simons coefficient on the different interactions and generalize 
the formalism to include various degrees of freedom arising in bilayer graphene which give rise to 
different topological phases. 
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I. INTRODUCTION 

Since its recent synthesis, graphene |l|, l2| has become an ideal playground to establish fundamental concepts in 
condensed matter physics [3( . It is by now well known that the low energy excitations of the neutral system can be 
described by a continuum model that includes two flavors of massless Dirac particles in two spatial dimensions [J] . 
The two flavors are associated to the two Fermi points that constitute the Fermi surface of the neutral material 
and, in the standard description, the electron spin is an extra degree of freedom usually disregarded as long as 
^ . the physics in the system does not explicitly depend on it. In the last years it has been realized that when the 
physics involves the spin degree of freedom, interesting things might happen and new states of matter can emerge 
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like the topological insulating state [5j-|8|]. This state relies on the existence of a nonzero spin orbit (SO) coupling 
which is believed to be very small in graphene [5J, [9( but which recent proposals show that it can be enhanced 
by various factors such as curvature [10, LUj j impurities [l2| or Coulomb interactions [13[ . Also the application 
of a perpendicular electric field in graphene can induce an extrinsic Rashba SO coupling whose value can also 
_JL ' be enhanced by interactions with the substrate [l4j . Apart from the SO processes it is known that a staggered 
r-t potential can open a gap in the spectrum, existing some experimental evidence about the gap opening process 
Q when graphene is deposited on top of SiC 15]. The way in which those terms compete among each other can lead 
O ■ to a rich phase diagram of insulating and gapless phases as discussed in [16j | for monolayer graphene. 

There are other ways to get a nontrivial insulating state in graphene without invoking the presence of a SO 

fvj , coupling. Recently it has been shown theoretically that a proper choice of strain pattern can lead to the appearance 

^ ■ of a gapped phase, with an energy spectrum similar to the sy_stem of Landau levels and a pair of counterpropagating 

edge states, carrying different value of the valley index [17J . Contrary to what happens when the SO coupling 

is considered, here we cannot establish that these counter-propagating states are protected by some symmetry 

(when the SO coupling is considered the edge states are protected against disorder by time reversal symmetry). 

; Nevertheless, the topological character of the counter-propagating states is maintained for distances smaller that 

C ' the mean free path in the system. Other situation where a topological insulating phase can appear without a SO 

^^ interaction is by deposition of ad-atoms on the surface of graphene. When these ad-atoms are placed only in one 

of the two sublattices, a gap is opened and again, a pair of counter-propagating states appear at the boundaries 

carrying opposite valley number [18j |. 

The goal of this work is to describe how time reversal invariant topological phases can be identified and analyzed 
in a unifying way using the effective action formalism, a technique which has been widely used in condensed matter 
systems [X9J] . We derive a way to compute the Chern-Simons coefficient associated to the different situations allowed 
C^ by the bilinear terms in the hamiltonian. We will also see how terms different from the usual spin orbit coupling 
can lead to a Chern-Simons term even in the absence of an external electromagnetic field. It is well known that 
in the graphene system a finite density of carriers can be induced by applying an external voltage. We will pay 
special attention to the situation in which there is a finite chemical potential in the samples. As we will see this 
changes some physical situations and introduce phases in which the system is metallic with a finite value of the 
Chern-Simons coefficient. 

Bilayer graphene [20( is a more promising material than the monolayer for potential applications due to the 
possibility of opening a gap by applying an external voltage [2l| . Its possible insulating phases are also richer and 
are being subject of great attention recently [22 - l25l |. We extend the formalism to the bilayer case and provide a 
table of the insulating phases arising from the competition of the various spin-orbit couplings discussed in [2J| . 

The article is structured as follows: In section [II] we describe the model for the monolayer graphene used in 
this work and present the formalism of the effective action to derive an expression to compute the Chern Simons 
coefficient. In section HTT1 the formalism is applied to study various insulating phases in the monolayer case with 



a finite chemical potential. Section [IVI is devoted to analyze various non-trivial topologically insulating phases in 
bilayer graphene. We conclude with a summary and discussion of the results in section IVI 

II. THE SPIN CHERN SIMONS COEFFICIENT FROM THE EFFECTIVE ACTION 

In this section we review some of the known low energy aspects of graphene in the context of the effective action 
formalism for completeness. For details, we refer the reader to reference [16[ and to the references mentioned 
below. 

A. Monolayer graphene: The model 

The low energy elementary excitations of neutral graphene around its Fermi points can be described by the 
following massless Dirac hamiltonian: 

"H (k) = i[<J y T z k x - <7 x l T k y ]l s = 7 1 fc E +j 2 ky. (1) 

The Pauli matrices labeled a, t and s represent pscudospin, valley and spin degrees of freedom respectively. The 
7 matrices are 8x8 matrices chosen to fulfill the condition {7^, 7"} — Ax 2g^ v , with g^" = diag(l, — 1, —1). With 
these conventions, the zeroth gamma matrix is defined as 7 = l s l v a z . Since we will not consider many body 
effects the Fermi velocity is set to one by a proper scaling of the fields and the parameters in the hamiltonian. We 
have also put K = 1. Note that the matrices chosen differ from the convention used in [5[ ("/km)- The two are 
related by 7^ = cr z j KM . In what follows we will group all possible bilinear couplings with an arbitrary matrix A 
and write the corresponding Lagrangean as 

c = <k 7 % - E A * A ^> ( 2 ) 

i 

where ijj — i/' + 7o- While each 2x2 Hamiltonian around a given Fermi point is not time reversal invariant ( T: 
H(Ki) = H*(K2) ) grouping the two valleys in a 4x4 notation makes the full Hamiltonian Ho invariant under 
both inversion symmetry (realized as an interchange of sub-lattices A and B often called parity V in the graphene 
context ) and time reversal T(interchange of K\ and K% ) . These discrete symmetries -together with translation 
invariance- protect the Fermi points and prevent the opening of a gap in the spectrum when the intrinsic spin of 
the electron is not taken into account |2q |. 

A staggered potential is a time reversal invariant contribution (although it is a Dirac mass term that breaks 
parity) that opens a gap in the spectrum. In our notation it will will be represented by a coupling constant m 
followed by the unit matrix 

A m = 1 S 1 T 1 CT . (3) 

The main contribution of [5j was to realize that when the spin degree of freedom is taken into account the 
intrinsic SO term acts like two copies of the Haldane mass 27], with opposite sign for each spin projection. The 
Haldane mass was one of the first attempts to get Landau levels with net zero magnetic field [28j. With our 
convention it is defined as Ah = 1 s t z 1<j and breaks T. The symmetry is restored by the inclusion of the two 
opposite spin contributions. In our convention the intrinsic SO coupling has the form: 

H so = A so A so = A so s z T z l a . (4) 

We will also consider in our work the Rashba term which is written as 

Hr = XrAr = ^Ri(syT z <j y + s x \ T a x ). (5) 

As mentioned before, the SO term opens a gap in the system and respects both time reversal and parity symmetry. 
The other time reversal invariant term that opens a gap is a spatially constant Kekule distortion which is off- 
diagonal in the valley and lattice indexes. Its contribution to the Chern-Simons term discussed in this work is the 
same as that of the staggered potential and we will not consider it here. 

B. The spin Chern-Simons term 

It is well known that in (2 + 1) space time dimensions, a single massive Dirac fermion coupled to a background 
U(l) electromagnetic (EM) field induces a topological Chern-Simons term in the effective action of the gauge field 



29, 30]. In condensed matter language and in connection with the physics of the quantum Hall effect [28, 31, 32]) 
this can be seen as follows: The action for the minimal model of a Dirac fermion ip coupled to an EM field described 
by the gauge potential A^ is given by 



S = d 



-fl . 



^7"(d„ + ieA^ + m^ ~ \f^F^ v 



(6) 



The first radiative correction to the vacuum polarization (related by the Kubo formula to the conductivity tensor) 
is given by 



n ^ (9)=le2 /(2^ Tr 



1 „ 1 



7 M -T7 — -r 



(7) 



7 p fc p — m "fp(k + q) p — in 
From the properties of the Dirac matrices it is known that the vacuum polarization tensor can be decomposed as 



IV((?) = ( 9lu , - ^ ) ir(<z 2 ) + *m W gTlV), (8) 



where the symmetric (antisymmetric) part is even (odd) under inversion and time reversal symmetries. The tensor 
structure of the odd term, which emerges from the trace of three 7 matrices, induces a Chern-Simons interaction 
in the effective action of the electromagnetic field of the form 



C / d 3 xe^A^A p , (9) 



and a transverse fermionic current: 



< jm >= Ce^ p d v A p , (10) 



where C is the Chern-Simons coefficient which is proportional to the sign of the fermionic mass. 

When two Dirac fermions related by time reversal symmetry are considered, each of them will contribute to C 
with an opposite value giving a null transverse current (|10p. Time reversal symmetry dictates that C^ = — Cj. and 
hence the total Chern number is C = Cf + C4. = 0. 

So in order to generate a Chern-Simons-like term in a time reversal invariant fashion we need to generate different 
masses (open a gap in the spectrum) for different members of the same pair of flavors. 

In graphene there are four flavors of massless Dirac fermions related by pairs (valley number) by time reversal 
symmetry. To describe an insulating state originated by some mass like term, we need to couple the Dirac 
fermions to a proper background field which can describe the difference between a standard insulating phase and 
a topologically insulating phase. Given a flavor, say spin, we must couple the fermions to an external background 
Vft field that sees the spin: 

JA^ (7 ^ _ A ty _ ^ A „ - H'-kWI (11) 

The matrix 75 has to be chosen appropriately to do the job and it should not be identified with the usual product 
of all the other matrices. Its role is to ensure that the traces in the calculation of the gauge effective action are 
nonzero (the details are explained in the next section). For instance, in the case of the spin current discussed in [5[ 
the matrix 75 is defined as 75 = s z l„l CT . In this way after integrating out the fermions and calculating the lowest 
order terms in the effective action (see next section) , eq. ([TU]) will be properly modified to give 

<Jf-Jj>=^=C.e'»''a,4, (12) 



t "i '- SV, 
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The spin Chern-Simons coefficient C s is defined as C s = CTj- — Cj. [5( and it will be related to the electron's Green's 
function in the next section. As mentioned before, time reversal symmetry dictates that the total Chern number 
is zero but we can still find a transverse current being topological in origin given by C s . This number classifies 
the quantum spin hall state. 

This approach can be generalized to other degrees of freedom such as the valley degree of freedom by properly 
modifying 75. We will discuss this issue in the section devoted to the bilayer case. 



C. The effective action 

We will describe in this section how to construct the spin Chern-Simons coefficient (or any other topological 
charge) on practical grounds when a generic mass like term A is considered in the Lagrangean @, leaving the 
particular cases for the following sections. A similar construction was done in [33[ in order to classify the topological 
insulators in three spatial dimensions and in 16] for monolayer graphene. Consider the action defined in dill) . 
The one-loop effective action for the gauge field is obtained by performing the functional gaussian integration of 
the fermionic variables. Since (|11[) is quadratic in the fermions it is straightforward to get: 

/d 3 k 
— Trln( 7 ^-A- e7 ^- 7 ^ 75 K). (13) 

It is important to note that in equation (|13[) the momentum variable fc M actually represents the derivative operator, 
so in order to deal with the logarithm in (fl"3)) we need to expand it in Taylor series of the potentials A^ and V v . 
For that, we will formally define the fermionic Green's function as 



G(k) = — -, (14) 



and we will write (fT5|) as: 



ri[A,V] = J ^53__Ti[(Go(*,A)Hey*^-*^57 , '7lsK)) B ], (15) 

where we have dropped out an irrelevant (infinite) constant that does not depend on the fields A or V. In order 
to extract the spin Chern-Simons term we can concentrate in the I — 2 term in this Taylor expansion. Moreover 
since we will consider k in (|15[) as the derivative operator, we will seek for the terms in the expansion that will 
contain derivatives of the fields A and V, particularly the term VdA, as it can be seen from (fT2")l . We can perform 
the calculation in momentum space where the I = 2 term reads: 



p l f rl A n rt A k 

T[=2 = 7 J WfWT Tr t G ( fc )^ G ( fc + ff)^^] VvA,>. (16) 

Expanding G{k + q) in series of q up to first order we can identify the corresponding Chern Simons term which 
in real space reads: 

S scs = tfd 3 xC^ p V^A p , (17) 



where C^ vp is defined as 



p 2 r jus 
Cr = -i\ J ^pTr [G(k, K)YG{k, A) 7 "G(fc, A)y 75 ] . (18) 



It is interesting to note that the the last expression can be written in the form of a Pontryagin index [3J, l35j : 

d 3 k 
2 J (27r) 



C'r r = ~i%r I j—^TriGdk.G^Gdk.G-'Gdk.G- 1 ^]. (19) 



In what follows we will use these definitions to calculate some particular cases including the effect of the chemical 
potential. 

III. TOPOLOGICAL INSULATING PHASES IN MONOLAYER GRAPHENE 

A. An illustrative case: Intrinsic spin-orbit coupling 

Although known in the literature for quite some time [5( it is illustrative to consider the case with only intrinsic 
spin orbit coupling as a starting point to construct more elaborate situations. In the simple case of the SO coupling 



with the chemical potential (/i) set to zero, it is trivial to check that A so commutes with all the 7 matrices and 
thus we can write G(k) as 



, 7 p fc p + A so A s 
A 2 



G{k) = l ^ L ^ L - (20) 



so 



Written G(k) in this way, by the trace of the 7 matrices it is apparent that only the term proportional to A|? will 
survive, being proportional to Sie^ pu as well. The integral in momentum can be evaluated easily, as long as we 
note that the integral is insensitive to the sign of A so . The result is [_|: 

C avp _ tivp±_ / "*- ^^so _2.ii/p ^so _ 2 e „;„»(*, \ (<y\ \ 

This is a widely known result. It can be proven that this number, sign(A so ) is the same Z2 topological invariant 
defined by Kane and Mele in the case of two dimensional topological insulators [6| . 

This result is modified by the presence of a finite chemical potential. When |/i| > A so , the chemical potential 
crosses one of the bands changing the position of the poles in (IT51) . Thus, we obtain: 

Cr = e^^r, (22) 

a result first obtained in [3_|. When \fj,\ < A so the two poles in (fT8|) are located always on different semi-planes 
(upper and lower) and so the integral gives the same result as in the case of zero chemical potential. With this 
result we learn how to proceed with more general cases and the effect that we should expect with the introduction 
of a chemical potential. When the chemical potential crosses the bands, the topological nature of the effective 
action is not destroyed, but the result is no longer quantized and depends explicitly on the position of the chemical 
potential. When the chemical potential lies inside the gap, the system is a topological insulator with a quantized 
spin Chern-Simons response. 

B. The competition between the intrinsic spin-orbit term and the staggered potential 

We will follow the method sketched in the previous section to study analytically the competition of the SO term 
with the staggered potential in the presence of a finite (and zero) chemical potential. The interest of this case 
resides on the fact that while both terms open a gap in the system, the nature of the insulating phase is different. 
It also introduces some technical complications whose resolution will be useful for the analysis of more complicated 
situations. The competition of the intrinsic spin-orbit term and the staggered potential was studied numerically 
in [5| for zero chemical potential. In what follows we will complete the discussion and give an analytical result 
valid for finite jjl. A similar analysis was also done in [l_| where the interplay of a Haldane mass and a staggered 
potential for zero chemical potential is analyzed. 

The appropriate mass-like term in this case is M = mA m + A so A so where A m and Aso are defined in ([_J and 
d_J respectively. 

We begin by the case where the chemical potential is zero. The first difficulty that arises is that we can no 
longer write the Green's function in the simple form (12U1) . In the simple case of the SO coupling all the fermion 
flavors had the same dispersion relation, w(k) = ±-\/k 2 + A 2 , and we were dealing with a multiflavor two band 
model. We only had an unique pole in the Green's function and the degeneracy of the system came from the matrix 
structure in the numerator of ([20)) . When considering the two terms, this degeneracy is partially lifted and we have 
a non-degenerated multiband model. It can be easily seen by diagonalizing the hamiltonian in the presence of the 
two terms. The four bands are given by the expressions u — ±-\/k 2 + (to + A so ) 2 and tu = ±-\/k 2 + (m — A so ) 2 . 
Effectively, our system is now made of two (doubly degenerated) two-band subsystems. However, although the 
denominator in (fT4)) is not any more proportional to the identity matrix, it is still a diagonal matrix that can be 
inverted. Using expression ([T5]) and following the recipe of the last section one can see that we have two copies of 
the SO problem with two different masses given by to ± A so . The result is accordingly: 



liup 



C £»p = e 2 -— [sign(A S0 - m) + sign( A so + to)] . (23) 

47T 

It is clear that the interplay between a staggered potential and a spin orbit coupling is such that the Chern-Simons 
coefficient is still quantized. When A so = we recover a topologically trivial insulator and when m = we recover 
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the result of the previous section as expected. The new feature compared to previous results is that when both 
are non zero (and for A| > \m\) there is still a topological response of the system. 

Lets turn now to the case of having a finite chemical potential, not discussed previously in the literature. Similarly 
to what happened in the intrinsic SO case, depending on the relative value of [i against A so ± m we will have 
different results. Without doing any extra work we can read the result from the considerations made in the case 
of the intrinsic spin-orbit coupling by changing the masses appropriately. There are 4 different cases: 



• 



\n\ < \A so - m\ and |/x| < A so + m\: 

Under these conditions the chemical potential lies inside both gaps gap and the response is still quantized: 

,p,vp 



• 



Cr P = e 2f — (sign(A so - m) + sign(A so + m)). (24) 

\n\ > \A so - m\ and |/i| > |A S0 + m\: 

In this case, the chemical potential crosses the bands for both masses and hence both integrals give a 
non-quantized result analogous to (J22I) : 

C , VP e 2^_ ( Ko-m A so + m \ 2 e^ / A^\ 

Cs ~ e 4. (, h + h )- e 2. UlJ' (25) 

It is interesting to note that this result docs not depend on the value of m. Finally, the two cases left can 
be written in a compact way: 

\n\ > \A so =p m\ and |/x| < |A S0 ± m\: 
These are the cases where one integral gives a topological contribution, as the chemical potential lies inside 
one of the gaps, but the other subspace gives a nonquantized contribution as the chemical potential crosses 
one of its bands. Hence, the result is: 

Cr = e 2 ^ (^P + si ^( A - ± m)) . (26) 

The ± signs indicate wether it is the conduction or the valence band which is crossed. 

From these simple analytic results one infers that there is a competition between the staggered potential and 
the intrinsic spin orbit, conditioned by the presence of the chemical potential. The inclusion of the finite chemical 
potential completes the known results for this case. 

C. Competition between a Rashba coupling and intrinsic spin-orbit coupling at finite chemical potential 

The competition of the intrinsic SO coupling and the Rashba contribution was already studied in the original 
reference [5[ for zero chemical potential. As the Rashba term by itself does not open a gap in the spectrum 
it was found that the topological insulating phase only exists for absolute values of A so bigger than |Ar|- The 
inclusion of a finite chemical potential makes the physical analysis more interesting. Unlike the other cases, this is 
more complex and a general analytical treatment with arbitrary /i is messy. A detailed calculation can be found 
in [36j; in what follows we aim to complete their analysis with several new analytical results together with an 
interpretation of the divergences that appear in this case. 

1. The case fi — — A so 

The band structure when these couplings are included is shown in figure [TJ As it is known in the presence of a 
Rashba term A/j the spin degree of freedom is no longer a good quantum number and we do not expect the spin 
hall conductivity to be quantized. After performing the trace in (fj"8)) . the integral left is given by: 

e 2 r dk 3 U\%(ko + A so )(ko(ko + A so )-k 2 x + kl)+A so (-k%+k 2 + A 2 so ) 2 ) 

Cr - Y / TgyfM*)^ ~. " " ~, '-, (27) 

J ( ' ((-k 2 + k 2 + AL) 2 - 4A|(fc + A so ) 2 ) 



where k — ±-\/k 2 + (Xr — A so ) 2 — XR — fi and k^ B = ±-\/k 2 + {Xr + A so ) 2 + Xr — fi are the poles. The simplest 
case (regarding its pole structure) is the case fi = —A so (see figd]), since for all values of the parameters there are 
two bands below (and two above) the chemical potential, therefore defining the same pole structure for all values 
of k. Since the chemical potential is placed at a special point of the band structure were the bands meet when the 
gap closes, we expect to have divergences which we interpret as the divergent contribution from the zero modes. 
There are eight possible cases depending on the sign of Xr, A so , A so + Xr and A so — Xr. Since the integral in 
(|2~T1) does not distinguish the sign of Xr we can safely assume A^ > 0. This reduces the problem to four cases: 

• A so > 0, A so + X R > and A so - X R > 0: 
In these conditions the gap remains open (A so > Xr) we get a finite result: 



cr p 



47T Xr 



-In 



Xr + A s 



Xr - A s 



(28) 



Note that when Xr/A so — > we can make use of the relation 

arctanh(x 



2 \l-x 
to recover the standard result with the sign function obtained for A so only. 

• A so < 0, A so + X R < and A so - Xr < 0: 

With an analogous analysis one arrives to the same expression but with a sign difference, which accounts for 
the change in sign of A so . 

• A so > 0, A so + X R > and A so - Xr < 0: 

In this case, one can see that the integral diverges. As explained above, we now have a semi metal with the 
chemical potential set at the touching point of the bands. The contribution of the k = particles turns the 
integral to be divergent. 

• A so < 0, A so + X R > and A so - X R < 0: 

Again, the integral is found to be divergent but with an opposite sign. 

We can summarize the above result in two compact expressions. For the case when the gap is opened, i.e. 
A so > Xr > and A so < —Xr < 0, one can write: 



Cr p = e 2 sign(A so 



e^P A s 
4tt Xr 



-hi 



A, 



A, 



A, 



(29) 



For the other case (A_r > A so > —Xr) where the gap is closed the integral diverges as C£ up —> sign(A so )oo 
as a result of the particular position of the chemical potential, where it encounters an infinite contribution from 
degenerate zero modes. 

2. The case fi — 

When fi = it is immediate to notice that two of the bands can cross the fermi energy depending on the value 
of the parameters. Hence, the chemical potential lies inside the gap or crosses one of the bands, depending on the 
value of the parameters (See figure[T|). Again the integral (|2"T)l which we have to evaluate has poles at the dispersion 
relation. An important difference with respect to other cases is that the poles have a non-trivial dependence on 
k. We consider some illustrative cases: 

• A so > 0, 2X R > A so > X R : 
This corresponds to a situation where the gap is open but the chemical potential is such that it intersects 
one of the bands, hence being the system metallic. The result is given by: 



Cr p = e 2 - =^ln 

4ttAr 



2(Aj 



(30) 



Interestingly, the case where A so > 0, Xr > A so has the same result. In this case the chemical potential 
intersects one of the bands but the gap is closed. These results suggest that when the chemical potential 
intersects one of the bands and hence the system becomes metallic, it does not matter whether a gap is open 
or not below the chemical potential as both situations give the result shown above. 





(a) 



(b) 



(c) 





(d) 



(e) 



FIG. 1: Band structure for different values of Xr and A ao . (a) With A_r = and A so /Oa gap opens. The bands are 
degenerate in spin, (b) With Xr 7^ and A so = the gap is closed but the spin degeneracy is lifted, (c) With Xr < A so a 
gap opens in the spectrum. (d) With Xr > A so the gap closes, (e) When Xr — A ao three of the bands touch which leads to 
a divergence in C^" p (see text) . The horizontal solid line indicates where the chemical potential would lie if it was fixed at 
jj, = — A ao . The band structure at this point generates a divergence when Xr > A so . 

• A so > 0, A so > 2X R : 
In this case, the chemical potential falls inside the bulk gap and so this situation is reminiscent of the 
situation with \x — — A so when the gap is open. The result is the one given in (|28[) . 

As expected, we recover the result for the previous section, but, since the chemical potential is not at a singular 
point we do not encounter divergences. The other cases left can be worked out in a similar fashion. Since the 
information they carry is redundant they will not be discussed further. 

We have thus completed the analysis made in [36] by introducing several novel analytical results together with 
an interpretation on the origin of the divergences. The competition of A SOi Xr and staggered potential at zero 
chemical potential was studied in [6] and we will not analyze it here. 



IV. BILAYER GRAPHENE 



A. The model 



One can extend the previous analysis to bilayer graphene where in addition to the valley (r), spin (s) and 
pseudospin (er) degrees of freedom we have the layer degree of freedom (fj,). The new degree of freedom enables a 
richer playground for topological phases. Remarkably, the effective action formalism allows to span these phases 
with little effort, as we shall discuss next. 

We will consider the bilayer graphene Lagrangean given by 



4 b) = <K 7 % - J2 A * A ^' 



(31) 



where the bilinear couplings will be described below and the 7 matrices are now constructed in the space of 
valley, pseudospin, spin, and layer. Following the notation of the monolayer they are chosen to be 70 = 1 T 1 S 1 M (7 Z , 
7i = irzlsl^Oy, 72 = -HAsl^x. 



The bilinear terms include the standard interlayer coupling given in our notation by 

A*A t = -t±l s l T i((J, y (T x - (J-xVy), 



(32) 



and an external gate potential which is known to open a gap in the bilayer system and which is given by 

Ay Ay = Vl s l T fi z a z . (33) 

In addition we shall consider two types of bilinear terms, the staggered potential and the various spin-orbit 
couplings given in [241 ]. The staggered potential is just the trivial extension of the monolayer. It is given by 

A^AW = ml T l fl l M l ff , 



where the superscript (6) denotes bilayer graphene. 

Following [24[ (see also [37|), one can construct the following spin orbit terms (in the notation of j5|): 

^it-A-it = ^lTzSzl^a-z + X 2 r z s z fj, z l a + X 3 (l T s x /j, z a v - r z s y iJ, z (x x ) + 

+ A4 (lrSxfJ'y&z — T z Syll x a z ) . 

The first term is the one corresponding to the (monolayer) intrinsic spin orbit coupling. The last three terms, 
which involve /j, are intrinsic to bilayer graphene where the ones proportional to A3,A4 are Rashba-like term. 
We will see that the effect of Ai and A2 is to open a non-trivial gap in the spectrum and so turning the 
system into a topological insulator with unprotected edge states. Rashba-like terms where studied numerically 
in (38( and shown to have the same effect as in the monolayer, and hence we will not study their effect in this work. 

The usual spin Chern-Simons term analogous to the one discussed in the monolayer case will be associated to 
the matrix 

75 = l T s 2 l M l CT . 

We will denote as before the value of the coefficient by C s . Its computation follows the lines described in sect. Ill CI 
A non zero value indicates that the system is a spin topological insulator and the edges support spin currents. 

As discussed in section [II] the present formalism allows us to describe other types of insulators by appropriately 
redefining 75 to resolve other degrees of freedom, which certifies the power of the method. In particular we 
will compute a "valley" Chcrn number to study the "valley" topological insulator character of the system. This 
insulator will be characterized by a matrix 



75,: 



T Z 1 S 1 M 1 C 



tailored to resolve the valley degree of freedom. The corresponding Chern-Simons coefficient will be named Ck ■ 
A non zero value of this coefficient implies a topological insulator with valley polarized edges. 

B. Summary of the results for bilayer graphene: 

We have computed the C s terms for bilayer graphene along the same lines described previously for the case of 
the monolayer. We summarize the results in Table I where the non zero couplings are listed (t± is always non 
zero and positive). As discussed above, we have also calculated the Ck Chern number which is the "valley" Chern 
number calculated with 75 replaced by 75^ in the trace (|I8p and get similar results. 

We omit an explicit derivation of the results since the integrals appearing are similar to the monolayer case, which 
we have detailed earlier and no new technical issue appears for the case of the bilayer. Nevertheless, the table also 
includes the technical detail of why the result vanishes in cases when it is zero, and the method of evaluation of 
the spatial integral if the result is finite. The chemical potential in all cases is set to zero, which in all cases falls 
inside the gap. With this set up, no band crosses the Fermi energy and the integrals are easily evaluated. 



Non zero couplings 






c s 


C K 


V/0 


(trace vanishes) 


^-sgn(V) (numerical) 


Ai/0 


^■sgn(Ai) (analytic) 


(trace vanishes) 


A 27 ^0 


^-sgn(A 2 ) (numerical) 


(trace vanishes) 


A3/O 


(vanishes by parity) 


(trace vanishes) 


A 47 ^0 


(trace vanishes) 


(trace vanishes) 


Ai / A 2 / 




2e 
It 


-sgn(Ai + A2) (numerical) 


(trace vanishes) 


Ai / V / 


~7 T 

2 77 


(sgn(Ai 


+ V) + sgn(Ai — V)) (numerical) 


f^(sgn(V + Ai) + sgn(V — Ai)) (numerical) 


A2/OV/O 


7.TT 


(sgn(A 2 


+ V) + sgn(Aa — V)) (numerical) 


f^r(sgn(V + A 2 ) + sgn(V - A 2 )) (numerical) 


ra/0 


(trace vanishes) 


^rsgn(m) (analytic) 



TABLE I: Classification of the two topological insulating phases of bilayer graphene considered in the text. 



These results indicate that bilayer graphene can be a "valley" topological insulator, as it has been recently 
suggested [39j by only applying a gate voltage V and also a topological insulator when Ai and A2 - which open 
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non trivial gaps - are present. However, note that the Chern number for bilayer graphene is two times bigger 
than for the case of the monolayer in all cases. This even Chern number implies that we are dealing with a 
topological insulator with unprotected edge states. The result that bilayer graphene is a topological insulator for 
Ai was already mentioned in [24j but nothing was said about A2 which also seems to make the bilayer system a 
topological insulator. 

The origin of a non-trivial gap given by A2 can be traced back to the matrix form of this coupling, together with 
the form of a t± coupling in bilayer graphene and the structure of 75 . All these together provide the complete set 
of matrices in the layer index (fi x , fJ. y ,n z , 1/*) so that the trace does not vanish and a non trivial spin Chern-Simons 
is generated, having however a "layer-like" origin. It is easy to check that when t± — 0, A2 does not contribute 
to the spin Chern-Simons term. In contrast, Ai in this situation still does contribute, indicating that the coupling 
between layers is of critical importance for this spin-orbit like coupling to have an effect. 

From the table of results it is clear that the combined perturbation Ai and A2 can enhance the topological 
response. The reason why the couplings Ai and A2 have the same effect can be understood in terms of the 
effective low energy hamiltonian. The equivalence can be traced back to a symmetry analysis of the effective 
model performed in ref. [33] . This analysis shows that the only term allowed by time reversal and discrete spatial 
symmetries is the Kane-Mele term associated to Ai. It means that although Ai and A2 have a different microscopic 
origin, both couplings will lead to the same term in the low energy approximation explaining why in table I the 
only effect of A2 is to renormalize Ai. 

The case with Aj j^ (i = 1,2) and V 7^ is analogous to the monolayer case where the staggered potential 
competes with the spin orbit coupling. As expected, the trivial coupling is V when calculating C s and Ai when 
calculating Ck ■ This indicates that the (valley) topological nature of the bilayer is affected by the presence of the 
gate potential V. 

As mentioned above, the Rashba like terms A3 and A4 go against the topological nature of bilayer graphene as 
it was discussed in (38| . Table 1 shows that bilayer graphene is a rich playground to understand the competition 
between different topological phases. The addition of a finite chemical potential will change this picture as it 
happens in the monolayer case. It is worth to mention that the computation of the different terms in table 
I is similar in many aspects to the computation of C s in the case of monolayer graphene in the presence of 
Xr + A so couplings. In the present case the existence of the perpendicular hopping term t± makes the low energy 
hamiltonian quadratic in the momentum operator and C s will have a form similar to the expression (|27p with the 
obvious modifications concerning the position of the Fermi level inside the gap. 

V. CONCLUSIONS AND OPEN QUESTIONS 

Wc have provided a unifying view to generate and compute different Chern-Simons terms that can appear in 
mono and bilayer graphene, and used the formalism to study various topological insulating phases in the system. 
Part of the results of this work were already obtained in the literature and others are distinct. In particular in 
the monolayer case we have studied the competition of the intrinsic spin orbit term with a staggered potential in 
the presence of finite chemical potential, and we have completed the analysis of [36[ by obtaining an analytical 
expression for the Chern-Simons term in the case of having an intrinsic SO and a Rashba term at finite chemical 
potential. We have also clarified the origin of the divergences encountered for the particular case of a chemical 
potential equal to the intrinsic SO value. For the bilayer graphene we have characterized the usual spin insulator 
and a different non-trivial topologically insulating phase similar to the one described in [39j, that can be thought 
of as a valley topological insulator, under the same formalism. The topological nature of the valley insulator was 
shown to be affected by the introduction of a gate voltage. We have also seen that the layer degree of freedom 
does not give rise to topological insulating phases with the physical couplings discussed in this work. Obviously 
the formal identity of the Hamiltonians with various couplings allows to play with possible couplings to trade spin 
by layer degrees of freedom. The major interest will lie in finding couplings that will be experimentally realizable 
and can lead to non-trivial applications as discussed in [39|, an issue out of the scope of the present work. 

The effective action formalism used in this work allows to compute non quantized values for the spin Chern 
numbers, in contrast to other methods, like the computation based on the Berry phase, which only gives the 
quantized parts. This non-quantized Chern number is still linked to a spin Hall effect which means that we still 
have counter propagating edge states whose spin polarization is given by this Chern number. A non quantized 
value means that the spin (or generally speaking the quantum number carried by the propagating edge state) is 
not strictly conserved. In spite of this, there have been some theoretical studies (in the case of valley Hall effect) 
claiming that although not protected by symmetries, these counter propagating edge currents are quite robust 
against the effect of disorder and the coherence length of the quantum number they carry can be significantly large 
to be used in some electronic or spintronic applications (39|. 
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